Limiting fragmentation of chemical potentials in heavy ion collisions 
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Thermal models have been used to successfully describe the hadron yields from heavy ion collisions 
at a variety of energies. For ^/sjvjv < 17 GeV this has usually been done using yields integrated over 
4-7T but at the higher energies available at RHIC, yields measured at central rapidity have been used. 
Recent BRAHMS data allows us to test whether thermal models can be generalized to describe the 
rapidity dependence of particle ratios. We have used the THERMUS package to fit BRAHMS data 
for the 5% most central Au+Au collisions for several rapidities at ^/sjvjv = 62 and 200 GeV. We have 
found a relationship between the strange and light quark chemical potentials, /is = 0.21 ± 0.01/ib. 
Using this relation we are able to describe the energy dependence of A, 3 and O ratios from other 
experiments. We also find that the chemical potentials are consistent with limiting fragmentation. 
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INTRODUCTION 

A wide variety of data from the four RHIC experi- 
ments is consistent with the idea that high energy heavy 
ion collisions produce an almost perfect partonic fluid |l| . 
The hot dense fluid expands in both the transverse and 
longitudinal directions until the quarks freeze-out into 
hadrons. In contrast to initial expectations these colli- 
sions produce particle distributions that are not boost in- 
variant along the beam axis but rather are Gaussian with 
widths that depend upon particle type, 0] ■ This suggests 
that the "chemistry" of the quarks may depend upon ra- 
pidity. Traditionally thermal analyses have assumed that 
there is only one source in heavy ion collisions and have 
used particle yields integrated over all phase space as in- 
put to grand canonical fits. At RHIC several groups have 
shown that such models can give an excellent descrip- 
tion of a large number of particle ratios at mid-rapidity 
d Q. The wide rapidity range of BRAHMS allows us 
to test if there is more than one thermal source. The 
THERMUS program is an implementation of the Grand 
Canonical Ensemble in the interactive C++ framework of 
ROOT [10. THERMUS assumes that at freeze-out all 
hadrons are in chemical equilibrium. Within this model 
the conservation of isospin, strangeness and baryon num- 
ber is controlled by chemical potentials and a freeze-out 
temperature T. The model assumes that the total net 
strangeness is zero and that the isospin of the system is 
set by the neutron/proton ratio of the colliding beams. 
While these two things are obviously true for 4-7T yields, 
they are assumptions in our analysis. 

BRAHMS has measured n , K ± ,p and p yields ver- 
sus rapidity at ^/JWn = 200 and 62 GeV, H Q 
(The 62 GeV data are still preliminary). From these 
data it is possible to construct five independent ratios, 
7r~/7r + , K~ / K + , K + /tt + ,p/p and p/n~ that are input 
to THERMUS. This is done for data at y=0,l,2 and 3. 
BRAHMS is a spectrometer experiment and so data for 
all particles of the same sign of electric charge is taken 
at the the same magnetic field setting. Therefore some 



interpolation is required to make ratios at the same ra- 
pidity. 



THE GRAND CANONICAL ENSEMBLE 

In a chemical analysis within the grand-canonical en- 
semble, the statistical-thermal model requires at least five 
parameters as input: volume V, the chemical freeze-out 
temperature T, baryon-, strangeness- and charge chemi- 
cal potentials ^b,HS and /zq respectively. In addition a 
strangeness saturation factor 75 is sometimes employed 
to allow for the possibility that the system will not have 
time to equilibrate strangeness. When considering parti- 
cle ratios the volume drops out. Since we only have five 
independent ratios we have to fix some of the parameters 
"by hand." For this analysis we fix take T=160 MeV and 
75 = 1. In addition /iq is constrained to reproduce the 
proton/neutron ratio for Au+Au collisions. This leaves 
five ratios to constrain two free parameters, /is and fj,B- 



FITS TO BRAHMS DATA 

Figure ^ shows fits of particle ratios measured in 
BRAHMS at ^/IWn =62 and 200GeV to THERMUS. 
The ratios are reasonably close but the values of x 2 /NDF 
are quite large, particularly at ^snn =62 GeV. This may 
indicate a failure of the model, an underestimation of the 
(preliminary) errors or the fact that we only have five par- 
ticles ratios to constrain the fit. For y=0 at ^/snn =62 
GeV, the fit is constrained by the requirement fig > 0. 
For this point we have redistributed the probability at 
negative /is to positive values using a Gaussian distribu- 
tion. This increases the error on /is for this point. 
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FIG. 1: (Top) The rapidity dependence of our tt~ /tt+ , K~ /K + ,K + /% + ,p/p and p/V" at v /5^" = 200 GeV (left) and 
preliminary 62 GeV (right). The data are shown by solid symbols and the Thermus fits by open symbols. (Bottom) the \ 2 P er 
degree of freedom. 
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FIG. 2: The chemical potential of strange quarks fis versus 
that of light quarks \ib for chemical fits to BRAHMS data 
using THERMUS. The line shows a least squares fit to these 
results. 
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FIG. 3: The chemical potentials versus y — ybeam- The upper 
line shows an exponential fit to the [is data. The bottom line 
is the same fit scaled by 0.21. 




FIG. 4: n/fi and S/3 ratios from STAR, NA57 and NA49 
versus predictions from THERMUS using fj,s = 0.21 • \ib @, 
[Il[lll . The 9, 62 and 200 GeV data are still preliminary. 



13 plots the results of these fits as a function of y — ybeam- 
Both chemical potentials rise as a function of y — ybeam 
and at a fixed rapidity decrease with energy. The upper 
line shows an exponential fit to the fis data. The data 
are reasonably consistent with the hypothesis of limit- 
ing fragmentation though the 62GeV data are somewhat 
above the curve while the 200GeV data are below it. The 
lower line is just 21% of the upper one. 



ENERGY DEPENDENCE OF HYPERON RATIOS 

Using this relationship THERMUS can predict the cor- 
relations between anti-hyperon/hyperon ratios and the 
p/p ratio since they are controlled, by the two fugaci- 
ties, e _MS / T and e~ MfJ / T . Note for these ratios that any 
strangeness saturation factor cancels out.) In Fig. 0] the 
f2/n, E/3, A/A ratios measured at y=0 for central PbPb 
and AuAu collisions at ^snn =9 - 200 GeV are plotted 
against the cor resp onding the p/p ratios from the same 
experiments |9l ll0Llll| . The curves show the correspond- 
ing predictions from THERMUS given n$ = 0.21 ■ (1b- 
The agreement between the prediction and data is good, 
although THERMUS is slightly below the fl/Q data. 
Note that to first order the uncertainty on the tempera- 
ture T does not affect the shape of these curves. 



CORRELATIONS & LIMITING 
FRAGMENTATION 

Figure |2 shows the correlation between fi$ and p,s- It 
is clear that there is a linear relationship between the 
two potentials. We have done a least squares fit and find 
lis = (0.21 ± 0.01) ■ fi B . with a x 2 /NDF = 34.3/7. Most 
of the x 2 is generated by the y=0 point at 62GeV. Figure 



CONCLUSIONS 

Using a Grand Canonical description of heavy ion col- 
lisions encoded in the THERMUS framework we can 
give a rough description of the rapidity dependence of 
hadron ratios from AuAu collisions at ^Jsnn =62 and 
200 GeV. Given that only five independent ratios are 
available there is a limit to how well we can test the ther- 
mal ansatz. Away from mid-rapidity we cannot prove 
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that T=160MeV or that there is zero net strangeness in 
each unit of rapidity. However if one accepts this ansatz 
several trends are evident. The chemical potentials of 
both light and strange quarks increases with rapidity 
and decrease with energy, implying that different regions 
of rapidity lose contact with each other before reaching 
chemical equilibrium. The chemical potentials seem to 
obey limiting fragmentation and rise exponentially with 
y — Ubeam- Finally the strange quark chemical poten- 
tial is proportional to the light quark chemical potential, 
[is = (0.21 ± 0.01) • [Lb- Using this relation we are able to 
describe the energy dependence of hyperon production. 
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